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ABSTRACT 

We present results from a 1 1 2 ks long look by XMM-Newton at the ultraluminous X-ray source 
(ULX) Holmberg II X-l, long thought to be the one of best candidates for the missing class 
of intermediate mass black holes (IMBHs). Our data comprises the first high quality XMM- 
Newton/RGS spectrum of an ULX, and an XMM-Newton/EPIC spectrum with unprecedented 
signal-to-noise. A detailed timing analysis shows that any variability on time-scales of minutes 
to hours is very weak (< few per cent fractional rms), though larger amplitude variations on 
much shorter time-scales could be hidden by photon counting statistics. This result suggests 
that if Ho II X-l harbours an IMBH, then we are observing this source in a highly unusual and 
atypical state when compared with the known variability behaviour of other accreting systems 
of large mass. Moreover unlike Galactic X-ray binaries, our spectral analysis indicates the 
possible presence of an optically-thick low temperature corona. Taken together our timing 
and spectral analysis suggests that the compact companion is most-likely a high luminosity 
analogue of black hole binary systems similar to GRS 1915+105, the Galactic microquasar, 
harbouring a compact object of mass no greater than 100 M . 
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1 INTRODUCTION 

Einstein and Rosat observations of nearby galaxies revealed a new 
population of point-like extra-nuclear X-ray sources with luminosi- 
ties far in excess of the Eddington limit for spherical accretion 
onto a 10 M Q black hole (L x > 10 39 erg s -1 , e.g. Fabbiano & 
Trinchieri 1987; Roberts & Warwick 2000), the so-called ultralu- 
minous X-ray sources (ULXs). In fact, the ~ 10 39 — 10 41 erg s _1 
X-ray luminosities of ULXs, if due to sub-Eddington rate spher- 
ical accretion onto a compact object, imply masses of ~ 10 2 - 
10 4 M0. ULXs may therefore represent the so-called "missing 
class" of intermediate-mass black holes (hereafter IMBHs), the 
purported primordial building blocks of super-massive black holes 
and galaxies (e.g. Ebisuzaki et al. 2001; Madau & Rees 2001). The 
best current evidence in favour of IMBHs underlying ULXs derives 
from recent XMM-Newton spectroscopy of several ULXs, where a 
"soft excess" in the spectra may be fit by a relatively low temper- 
ature (~ 150 eV) accretion disc model (e.g. Miller et al. 2003; 
Miller, Fabian & Miller 2004). This is considerably cooler than 
found in stellar-mass black hole X-ray binaries (XRBs), and points 
to a larger, cooler disc surrounding a more massive (~ 1O 3 M0) 
black hole. However, X-ray spectral results based on continuum 
modelling may be open to some ambiguity, so these detections 
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are not necessarily conclusive in all cases (Roberts et al. 2005). 
Other possibilities involving stellar-mass black holes have been 
mooted to explain the extreme X-ray luminosities of ULXs, no- 
tably anisotropic radiation and/or relativistic beaming (King et al. 
2001 ; Kording et al. 2002) or truly super-Eddington accretion discs 
(Begelman 2002; Ebisawa et al. 2003). 

A more direct determination of the mass of the compact ac- 
cretor in ULXs is the key to further progress in this field. Perhaps 
the best method for determining the mass of the compact object 
in ULXs from X-ray data is through a comparison of their Power 
Spectral Density (PSD) with those of XRBs and Active Galactic 
Nuclei (AGN). The PSD of all known accreting systems are re- 
markably similar (Edelson & Nandra 1999, Uttley et al. 2002). In 
AGN, PSD shapes corresponding to both low 'hard' accretion state 
systems (NGC 3783, Markowitz et al. 2003) and high 'soft' accre- 
tion state systems (NGC 4051, M c Hardy et al. 2004) previously 
seen in XRBs are now seen. Significantly, the location of the char- 
acteristic PSD break-frequency scales inversely with the mass of 
the black-hole (Uttley et al. 2002, Markowitz et al. 2003, Vaughan 
et al. 2005) for objects with masses ranging from 10 M©for a typ- 
ical stellar mass black hole (eg. Cyg X-l, tVcak ~1 Hz) to 10 - 
1O 8 M0 for the central supermassive black hole in nearby Seyfert 1 
galaxies. For the latter, black hole masses have been independently 
verified using both the stellar velocity dispersion and reverberation 
mapping techniques (Gebhardt et al. 2000, Ferrarese et al. 2001). 
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This method has been applied to ULX data on at least two previous 
occasions (Cropper et al. 2004, Soria et al. 2004); though putative 
breaks were detected, the data were insufficient to place strong con- 
straints on the black hole mass in either case. 

In order to progress, one requires the best possible X-ray data 
from an ULX. Here we present an X-ray spectral and timing analy- 
sis for an unprecedented long look with XMM-Newton at the ULX 
Holmberg II X-l (hereafter Ho II X-l). This is one of the most lu- 
minous (Lx ~ 1 - 2 x 10 40 erg s" 1 ) nearby ULXs, located in 
the d ~ 3.2 Mpc dwarf star-forming galaxy Holmberg II. Of the 
87 candidate ULXs catalogued by Colbert & Ptak (2002), Ho II X- 
1 had the highest observed ROSAT count rate making it an ideal 
target for PSD analysis. It is a very well-studied object; ROSAT 
and ASCA observations (Zezas, Georgantopoulos & Ward 1999; 
Miyaji, Lehmann & Hasinger 2001) revealed a relatively soft X- 
ray spectrum, and variability consistent with an accreting black 
hole. The detection of a soft excess in the X-ray spectrum, that 
can be modelled as the accretion disc around a ~ 1O 3 M0IMBH, 
was confirmed in three short (~ 10 ks) XMM-Newton observations 
(Dewangen et al. 2004). This ULX is also embedded in a spectacu- 
lar nebula, visible in both optical emission lines (Pakull & Mirioni 
2003, Kaaret, Ward & Zezas 2004, Lehmann et al. 2005) and radio 
continuum (Miller, Mushotzky & Neff 2005). Crucially, the flux 
in the high-excitation Hell A4686 line around the position of the 
ULX implies a lower bound on its isotropic X-ray luminosity of 
4 - 6 x 10 39 erg s" 1 (Kaaret, Ward & Zezas 2004). Hence Ho II 
X-l is one of the strongest candidate IMBHs in the local universe. 



2 A LONG LOOK AT HO II X- 1 WITH XMM-NEWTON 

Ho II X-l was scheduled for a 112-ks XMM-Newton observation 
on 2004 April 15 (observation ID 0200470101), though the onset 
of intense background flaring in the later stages of the observation 
caused it to be truncated to no more than ~ 80 — 95 ks exposure 
per instrument. The pn was operated in full-field mode, whilst both 
MOS detectors utilised the large window mode, and the medium 
filter was used for all EPIC detectors. Both the RGS and OM were 
in standard modes, in the latter case monitoring Ho II with a series 
of UVW1 observations (though we defer discussion of this data to 
a later work). 

Data for the analysis were extracted from the pipeline product 
event lists using the appropriate XMM-Newton SAS version 6.0.0 
tools. The pn data were filtered to leave only events with flag=0 
and pattern < 4; MOS data were filtered for pattern ^ 12 and the 
#XMMEA_EM flag. We quantified the background flare contami- 
nation by examining a 10 - 15 keV full-field light-curve from the 
pn detector. This showed the first ~ 45 ks of the observation to be 
relatively clean, with only a few minor flares. However, subsequent 
to that strong background flaring dominated the light-curve. In fact, 
a good time intervals (GTI) file selected from this light-curve for 
quiescent background times (count rate < 0.75 count s" 1 ) leaves 
only 38. 1 ks clean data for the whole observation. 

2.1 Spectral analysis 

The GTI file was used to produce the cleanest possible imaging 
and spectral products for the dataset. Ho II X-l is clearly detected 
in the data, with a combined EPIC (0.3 - 10 keV) count rate of 
~ 4 count s . EPIC spectra and response matrices were gener- 
ated using the ESPECGET tool. MOS and pn source spectra were all 
extracted from an 11 -pixel (= 44 arcsec) radius circular aperture 



over the position of Ho II X- 1 1 . A background spectrum for the pn 
was extracted from a slightly larger circular aperture on an adja- 
cent chip, at a similar distance from the chip readout nodes as Ho 
II X-l. MOS background data were taken from two small circu- 
lar apertures at the bottom-right and -left of each large window. 
RGS source spectra, background spectra, and response matrices 
were generated using the standard SAS pipeline script, RGSPROC. 
All spectral data were grouped to a minimum of 20 counts per bin, 
in order to provide adequate statistics for spectral fitting using x' 2 
minimization within XSPEC version 11.3.1. 



2.1.1 RGS observations 

The data extraction produced a total of 7088 source counts for both 
RGS modules combined, integrated over the whole (flare-filtered) 
exposure. The resulting spectrum and residuals, compared to the 
best-fit absorbed power-law model, are shown in FigureQ We note 
that this is the first reasonable quality, high resolution X-ray spec- 
trum of an ULX. A fit statistic of \ 2 = 430.4 for 367 degrees of 
freedom (dof) is obtained for the absorbed power-law model, where 
the absorption is modelled by a solar-abundance WABS model 
(based on the absorption cross-sections of Morrison & McCam- 
mon 1983), with a rejection probability (P re j) of 99%. The photon 
index over the 0.3 - 2.0 keV band was F — 2.6 ± 0.2, with an 
absorption column of JVh = 1.9 x 10 21 cm -2 , significantly above 
the known Galactic column density towards Holmberg II of Nn = 
3.41 x 10 20 citT 2 (Dickey & Lockman 1990). The measured flux 
over the 0.3 - 2.0 keV RGS band is 3.7 x 10~ 12 ergs cm" 2 s" 1 
corresponding to an unabsorbed luminosity of 1.2 x 10 40 ergs s _1 
at 3.2 Mpc. 

The residuals to the absorbed power-law fit show a clear ex- 
cess of counts between 0.5 - 0.6 keV, blue- wards of the neutral O 
edge. This excess can be fitted with a single, broadened emission 
line centered at 569 ± 9 eV, with an equivalent width of 28 ± 12 eV 
and a velocity width of a = 17 ± 5 eV (8900 kms" 1 ). The fit 
statistic then improves to x 2 /dof = 408.0/364. The line energy 
is consistent with emission from a blend of lines from the He-like 

VII triplet (at 561 eV - 574 eV). We therefore refitted this excess 
with three narrow, unresolved lines corresponding to the forbid- 
den, intercombination and resonance lines from the O VII triplet. 
Two lines are detected, the forbidden line measured at 563 ± 2 eV 
and the resonance line at 577 ± 2 eV, both with equivalent widths 
of ~ 7 eV, x 2 = 410.5 for 363 dof. The flux in either line is 
9 ± 5 x 10~ 5 photons cm" 2 s _1 corresponding to a line luminos- 
ity of ~ 1 x 10 38 erg s" 1 . 

One plausible explanation for the origin of the line emission is 
from an optically-thin thermal plasma, as has previously been sug- 
gested for this source (Dewangan et al. 2004). If we fit the MEKAL 
model in XSPEC to the spectrum, in addition to the underlying ab- 
sorbed power-law continuum, then we can place a limit on the con- 
tribution of such a plasma to the soft X-ray emission in HoII X-l. 
For solar abundances, it is found that the 90% confidence limit on 
the plasma luminosity is < 3.6 x 10 38 erg s _1 , consistent with the 
value found by Dewangan et al. (2004), whilst the plasma temper- 
ature is constrained at kT = 175 ± 25 eV. Thus an optically thin 
plasma is not likely to contribute more than 3% of the total soft 
X-ray flux in Ho II X-l. 

1 We limit ourselves to this relatively small aperture to avoid extracting 
data over pn chip gaps, though we note that > 87% of the source flux lies 
within this radius. 
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Table 1. XMM-Newton EPIC spectral fits. 
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Figure 1. Soft X-ray spectra of Holmberg II X-l obtained with the Reflec- 
tion Grating Spectrometer (RGS) on-board XMM-Newton. RGS 1 data are 
shown in black, RGS 2 in grey. The model fitted to the data (solid line) is 
an absorbed power-law, with solar abundances, as described in the text. The 
residuals to this model (lower panel) show an excess of counts between 0.5 
and 0.6 keV. 



Alternatively, the excess of counts above 0.5 keV may arise 
from incorrectly modeling the neutral O edge in the absorber 
model. To test this assumption, we refitted the absorbed power- 
law model using the TBVARABS absorption model, based on 
the improved abundances and cross-sections tabulated by Wilms, 
Allen & McCray (2000). We find that the O abundance is now 
0.56 ± 0.29 with respect to solar, while the column density (in 
excess of the Galactic value of 3.4 x 10 20 cm~ 2 ) is now Nh = 
1.6 ± 0.3 x 10 21 cm~ 2 . The fit statistic is improved significantly, 
with x 2 /dof = 384.8/365, i.e. a rejection probability of ~ 77%, 
so formally acceptable. The decrease in \ 2 for this fit is due to both 
the lower O abundance value relative to solar and the lower abso- 
lute abundance of O tabulated by Wilms, Allen & McCray (2000). 
The underlying photon index is then slightly flatter (T — 2.4±0.2), 
whilst the absorption corrected luminosity in the 0.3 - 2.0 keV band 
is 1.0 x 10 40 erg s _1 . In this scenario no O VII emission is formally 
required, as the depth of the O I edge is reduced, resulting in little 
excess emission between 0.5 - 0.6 keV. That Ho II X-l lies in a 
low-metallicity environment should perhaps come as no surprise; 
the O/H ratio for the Holmberg II galaxy is ~ 10% of the solar 
abundance (Hidalgo-Gamez, Sanchez-Sacedo & Olofsson 2003) 2 . 
Indeed, Soria et al. (2005, and references therein) suggest that low 
metallicity environments could favour the formation of larger black 
holes (w 50 M Q ), and hence preferentially host ULXs. Sub-solar 
abundance absorption along our line-of-sight to Ho II X-l there- 
fore seems the most likely explanation for the apparent 0.5 - 0.6 
keV excess in the RGS spectrum. 



2.1.2 EPIC observations 

EPIC spectral data were fit over the 0.3 - 10 keV band, though 
within this we only fit the MOS over 0.3 - 6 keV (above this 
the data were noise-dominated), and the pn over 0.7 - 10 keV 



2 The somewhat higher metallicity suggested by our RGS data in the local 
environment of Ho II X- 1 may be a consequence of the local star formation 
and/or the progenitor star for the compact object in the ULX. 
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The models used (subject to the tabulated absorption models) were, in 
XSPEC syntax: (1) PO; (2) DISKBB + PO; (3) DISKPN + COMPTT. 
We quote errors at the 90% level for one interesting parameter. The 
parameters are: a Absorption column external to our Galaxy in units of 
10 21 atom cm -2 . b Power-law continuum photon index. c Inner-disc 
temperature (keV). d Disc model normalisation. e Fraction of observed 
0.3 - 10 keV flux in the accretion disc component, f Maximum 
temperature of the disc (keV). 9 Coronal temperature (keV). h Coronal 
optical depth. 



(known calibration inaccuracies between 0.5 - 0.7 keV are par- 
ticularly prominent in this spectrally-soft dataset). We allowed for 
residual calibration uncertainties between the three detectors in our 
fits by introducing a constant term into our spectral fitting; in prac- 
tice the variation between the detectors was < 6%. Each model we 
fit to the data was subject to a fixed Galactic absorption column 
of 3.41 x 10 20 cm -2 , modelled using a solar-abundance XSPEC 
WABS cold absorption model. 

We began, as with the RGS analysis, by fitting an absorbed 
power-law continuum model, using a second WABS component to 
fit additional line-of-sight absorption to Ho II X-l (which is pre- 
sumably intrinsic to either the ULX itself, or its host regions in 
Holmberg II). This fit was rejected at high significance (see Table 
1), though its best fit parameterisation was similar to that of the 
narrower-band RGS result (T ~ 2.8, jV H ~ 1.6 x 10 21 cm" 2 ). 
However, we also attempted a second fit, this time using the 
TBVARABS model, with the actual abundance of the absorbing 
medium set to 0.56 times solar as suggested by the RGS O I edge 
feature. This produced a vastly improved fit (A^ 2 = 109 for the 
same dof). This improvement was predominantly down to the use 
of the TBVARABS model itself, which includes (amongst other 
improvements) a much lower absolute O abundance than WABS. 
We tabulate the results of these two fits, plus a TBVARABS fit 
set to solar abundance, in Table 1. However, even though the fit is 
greatly improved by the use of TBVARABS, a simple power-law 
continuum fit is still significantly rejected. 
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We next attempted two-component models consisting of a 
power-law continuum plus an accretion disc model, similar to pre- 
vious work on ULXs, and the empirical model for X-ray emission 
from a black hole X-ray binary. As with previous work, we utilised 
the XSPEC DISKBB model, based on the multi-colour disc black- 
body approximation (e.g. Mitsuda et al. 1984). This model provides 
an improvement to the fit; however, as Table 1 shows the magnitude 
of this improvement is highly dependent upon the choice of absorp- 
tion model, with the solar abundance WABS fits showing the largest 
improvement (A^ 2 = 158 for 2 extra dof), and the best actual fit, 
whilst the most physical absorption model (sub-solar TBVARABS) 
shows the smallest improvement (A\ 2 = 29 for 2 extra dof) and 
worst fit. In all cases the DISKBB temperature was relatively cool 
(~ 0.2 — 0.27 keV), albeit hotter than found in the previous obser- 
vations of Ho II X-l by Dewangen et al. (2004), and the power-law 
continuum remains steep (T ~ 2.6) 3 . 

The DISKBB + power-law model is only a first-order ap- 
proximation of the physics of an accretion disc plus corona 
system. Hence we next tried a more physically-realistic, self- 
consistent model, namely the DISKPN + COMPTT combination 
(see Gierlihski, Machiolek-Niedzwicki & Ebisawa 2001; Titarchuk 
1994 for descriptions of these models) that describes an accre- 
tion disc modified for e.g. relativistic effects, and a Comptonised 
corona. We link the two models by setting the seed photon temper- 
ature of the corona, To, to the maximum temperature of the disc 
T max , and we set the inner radius of the disc to the innermost 
stable circular orbit (= 6R g ). This model provides a statistically- 
acceptable solution to the spectrum, regardless of the absorption 
model used (as the modelled absorption is relatively low, there is 
little difference between the models in this case). Whilst the disc 
again appears cool (Tmax ~ 0.2 keV), the best-fit model shows a 
remarkable corona, in that it appears both quite cool (kT ~ 2.7 
keV) and optically thick (r ~ 4.4). We show the best fit spectrum, 
with the contributions from the two components decomposed, in 
Fig. [2] and tabulate the interesting parameters in Table 1. These 
models provide an observed 0.3 - 10 keV flux measurement of 
~ 6.7 x 10~ 12 erg cm -2 s _1 , which converts to an intrinsic lu- 
minosity of 1.0 x 10 40 erg s" 1 for a distance of 3.2 Mpc (this 
intrinsic luminosity is similar to the RGS estimate - taken from a 
much narrower band - as the absorption correction is considerably 
smaller in our more complex EPIC modelling). 

We have also investigated whether an optically-thin corona 
also provides a reasonable solution to the spectral data, as one 
might expect given the DISKBB + power-law fits to the data. This 
was done by attempting two new fits with the coronal tempera- 
ture set to 50 and 100 keV respectively, indicative of the range 
seen in Galactic black hole binaries, in the DISKPN + COMPTT 
models. These fits did indeed provide reasonable solutions to the 
data, with a best-fitting r < 0.25 (the other parameters did not 
change substantially from the values in Table 3 for the optically- 
thick solution) and P re _, ~ 95%. However, they were still worse 
than the optically-thick solutions by A^ 2 ~ 12 — 16 for one ex- 
tra dof, which equates to a significance of 3 — 4a according to the 
F-statistic. Therefore we regard the presence of an optically-thick 



3 We note that we also attempted the "non-standard" form of the DISKBB 
+ power-law model, as described by Stobbart et al. (2004) and Roberts et 
al. (2005). This model was only marginally less successful than the stan- 
dard usage, by Ax 2 ~ 15 for TBVARABS models, with best fits showing 
kTi n ~ 2 keV for the accretion disc, and r ~ 3 for the power-law contin- 
uum. 
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Figure 2. EPIC spectral data for Ho II X- 1 , plotted with the best-fit spectral 
model (3, with sub-solar abundance) and Ax-residuals. The pn data are 
shown in black, with MOS1 and MOS2 in grey. The contributions of the 
DISKPN and COMPTT components are shown by the dotted and solid lines 
respectively. 

corona as far more likely than an optically-thin one, though we ob- 
viously cannot exclude the latter. 

Finally, as Fig.[2]demonstrates, the data show no obvious signs 
of emission lines including Fe K. We place a limit on the possible 
contribution of a narrow, neutral (i.e. 6.4 keV) Fe K line by adding 
a Gaussian line model to the best-fit spectra; from this we derive a 
90% upper-limit of 25 eV on the equivalent width of such a feature. 

2.2 Temporal analysis 

Temporal data for Ho II X-l were extracted from the EPIC de- 
tectors in the same circular apertures as used in the spectral anal- 
ysis. Background data were extracted from same-sized apertures 
on the same CCD chip as the source. Given the relative spectral 
softness of the ULX, and the comparatively hard X-ray nature of 
background flaring detected by XMM-Newton, we maximised our 
temporal signal-to-noise by selecting only data in the 0.3 - 6 keV 
range. The use of relatively small apertures, in combination with 
excluding hard X-ray data, ensured that the source signal per unit 
area dominated above even the strongest background flaring. We 
demonstrate this for the pn data in Fig. [3] which shows the full 
light-curve for this observation of Ho II X-l binned to 100 s reso- 
lution. 

Although Ho II X-l shows substantial variability on time- 
scales of months-years (e.g. Dewangan et al. 2004), the XMM- 
Newton light curve shows very little variability on time-scales of 
less than a day (barring a very weak monotonic trend throughout 
the observation). To quantify and constrain the variability more pre- 
cisely, we extracted a 0.3-6.0 keV EPIC-pn light curve at the high- 
est possible resolution (73.4 ms for full-frame mode), correspond- 
ing to the first 48.2 ks of data, which is free of telemetry dropouts 
and significant background flaring. We split the light curve into 
five equal segments and measured the PSD of each segment (with 
power determined in units of (rms/mean) 2 ), then merged and av- 
eraged the resulting PSDs by binning in frequency, so that each 
bin corresponds to a factor of 1.1 increase in frequency (with the 
proviso that each bin contains at least 20 individual PSD measure- 
ments). The errors on the averaged powers were calculated from the 
rms scatter on individual PSD measurements in each frequency bin. 
The resulting PSD was found to be completely flat, consistent with 
variations being due to photon counting statistics only. To constrain 
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Figure 3. EPIC pn 0.3 - 6 keV light-curve for the long-look at Ho II X-l. 
The data is binned to 100-s intervals, and background subtracted. We also 
show the relative intensity of the background emission in the same band, for 
the same aperture. The gaps in the data are caused by the telemetry buffers 
filling up due to the high background flaring rate. For clarity we only show 
data bins where the we have 75-s or more of exposure. 



the strength of any intrinsic source variability, we used XSPEC to fit 
the PSD. We obtain an excellent fit (\ 2 = 74.3 for 88 degrees of 
freedom) for a constant power of 0.793 ± 0.002 Hz -1 (errors are 
90% confidence). The expected power level due to Poisson noise 
for the observed count rate is 0.8 1 Hz~ 1 . The fact that the observed 
constant power is smaller than expected may be related to weak 
instrumental effects (eg. pileup, deadtime). 

The fit residuals relative to the best-fitting constant are plotted 
in Figure |4] No systematic features are apparent. At first glance, 
the data might suggest that the fastest intrinsic variability of the 
source occurs on time-scales longer than observed (>day), thus 
supporting the idea that the source is an intermediate mass black 
hole. However, AGN with black hole masses ~ 10 e M show sig- 
nificant variability on time-scales of minutes-hours (e.g. M c Hardy 
et al. 2004), and we would therefore expect the same if the source is 
an intermediate mass black hole, unless it occupies a state with un- 
usually low variability power. To demonstrate this fact more quan- 
titatively, we attempted to fit the PSD with a constant plus a bro- 
ken power-law of low-frequency slope —1, and slope —2 above 
the break, similar to the PSD shape observed in several AGN and 
Cyg X-l in the high/soft state (e.g. NTHardy et al. 2004, 2005; 
Uttley & M c Hardy 2005). In units of powerxfrequency (vP(v)), 
a PSD slope of — 1 corresponds to a constant value Ci// (which 
when multiplied by natural log of 10 gives the power per decade in 
frequency). In the AGN observed to date, and in Cyg X-l and other 
black hole XRBs, this constant value is ~ 0.005 - 0.03. 

For a break frequency of 10 -4 Hz, which is even lower than 
those observed in several AGN (and corresponds to a black hole 
mass of ~ 10 6 Mq, assuming linear scaling of characteristic time- 
scales from the ~ 10 Hz break seen in Cyg X-l), the 99 per cent 
upper limit to Ci// is 2.8 x 10~ 4 . For a PSD break at 1 Hz or 
higher, which would correspond to black hole masses of 100 Mq 
or less, the limit is Ci// < 1.4 x 10~ 4 , which corresponds to an in- 
trinsic fractional rms over the entire observed PSD of less than 1 per 
cent. In some black hole XRBs in the high/soft state, the amplitude 
of variability is diluted significantly by the presence of strong con- 
stant blackbody emission from the disc which dominates the X-ray 
spectrum and dilutes the variability in the power-law component. 
However, our spectral fits show that the X-ray emission we observe 
is not dominated by the putative disc component (that accounts for 
only ~ 20% of the observed counts in the 0.3 - 6 keV pn data), 
hence we conclude that the source cannot be accreting in a state 
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Figure 4. Residuals from the observed PSD fitted with a constant. 



with timing properties similar to the high/soft state of Cyg X- 1 or 
AGN with good PSD measurements. 

We next consider the possibility that the PSD is 'band-limited' 
(e.g. see van der Klis 2005), that is, the source does vary signifi- 
cantly but that variability is confined to a limited frequency range 
(perhaps a decade or two, or even smaller if the variability is dom- 
inated by strong QPOs). Such band-limited PSDs are observed 
in black hole XRBs in the low/hard state (where relatively hard 
(r < 2) power-law emission dominates the X-ray spectrum) and 
also the very-high state, where the spectrum shows strong, steep 
power-law emission, similar to that observed here (see McClintock 
& Remillard 2003 for more detailed descriptions of these states). 
The high-frequency cut-offs in these PSDs typically occur around 
1-10 Hz, and simple mass-scaling arguments then imply that the 
peaks in band-limited power should be sampled by our observed 
frequency range, for intermediate mass or stellar mass black holes. 
In that case, we may not detect these peaks due to lack of sensitiv- 
ity, which decreases at high-frequencies due to the Poisson noise 
contribution to the PSD 4 . 

To determine upper limits on the amplitude of band-limited 
noise which may be hidden in the PSD, we fitted the PSD with 
a constant plus Lorentzian model. The Lorentzian represents the 
band-limited contribution: it's coherence or Q-value (the centroid 
frequency divided by the full width at half-maximum) can be 
changed to represent different observed types of band-limited noise 
or QPOs. For the incoherent band-limited noise typically observed 
in black hole XRBs in either state, we assumed Q — 0.5 (e.g. see 
Pottschmidt et al. 2003). For the strong, sharp QPOs which are also 
often observed in very-high states (and which would be easier to 
detect than band-limited noise of the same fractional rms), we as- 
sumed Q — 5 (e.g. see Remillard et al. 2002, van der Klis 2005). 
For these two types of Lorentzian, we fitted the PSD to obtain the 
90 and 99 per cent upper limits on variability amplitude (in terms of 
the square of fractional rms contributed by the Lorentzian), which 
is permitted as a function of Lorentzian centroid frequency. The re- 
sulting contour plots are shown in Fig. |5ja,b). The limits we can 
place on the presence of significant variability from band-limited 
noise are very stringent. For example, the typical observed frac- 
tional rms in a single Lorentzian is ~ 10 per cent, i.e. fractional 
rms-squared ~ 0.01, which is ruled out at 99 per cent confidence 

4 Note that for a band-limited component with constant rms and constant 
width in decades, the peak power also remains constant in uP(u) units, 
but the noise level increases linearly with frequency, hence the reduction in 
sensitivity at high frequencies. 
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Figure 5. Contour plots showing upper limits on the amplitudes of broad (Q 
PSD (see text for details) 

for frequencies less than ~ 0.05 Hz and < 1 Hz for Q = 0.5 and 
Q — 5 respectively. 

Assuming linear scaling of time-scales with black hole mass, 
these limits on the allowed frequencies of the expected band- 
limited noise suggest that the black hole mass should be of order 
100 M Q or less, assuming that X-ray timing behaviour is similar 
for all black holes, regardless of black hole mass. In fact the lack 
of observed variability is consistent with the PSD being similar to 
that observed in certain states of the microquasar GRS 1915+105, 
which, although famous for showing highly variable and complex 
light curves, also shows states (the \ c l ass of behaviours, see Bel- 
loni et al. 2000) where the long-term X-ray variability (on time- 
scales up to months) can be remarkably constant, with all the vari- 
ability power being concentrated at high frequencies (~ 1 Hz or 
higher). 



3 DISCUSSION 

Our analysis of the RGS data for Ho II X- 1 shows that it is likely 
that the X-ray emission we observe from the ULX is subject to ab- 
sorption by a sub-solar abundance medium, almost certainly within 
the Holmberg II galaxy, and possibly in the immediate environment 
of the ULX itself. Interestingly, this result has implications for pa- 
rameterisation of the EPIC X-ray spectrum. In particular, the mis- 
modelling of the neutral O edge as an solar-abundance absorber 
leads to an over-estimation of the magnitude of the soft excess 
in the empirical models (power-law continuum and power-law + 
DISKBB models) of the data. The main implication of this is that it 
can lead to an over-estimate of the mass of the putative IMBH de- 
rived from the DISKBB parameters. Using equations (5) and (8) of 
Makishima et al. (2000), and the best-fit parameters from the em- 
pirical power-law + DISKBB fits, the IMBH mass is estimated as 
624 [cos(i)]~ 5 a -1 M Q for the solar-abundance WABS absorption 
modeU. When the sub-solar TBVARABS model is used instead, 
the black hole mass reduces to 206 [cos(i)]~ 3 oT x M Q , down to 
~ 33% of the previous estimate. This explicitly demonstrates the 
sensitivity of the estimated IMBH masses to the assumed absorp- 
tion model, and may imply that previous masses in the literature 
derived by this technique are somewhat over-estimated. 

5 i is the inclination of the accretion disc to the line-of-sight, and a relates 
to the spin of the black hole (a = 1 in the Schwarzschild metric, and a < 1 
for Kerr black holes). 
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0.5, left) and narrow (Q = 5, right) Lorentzian components in the observed 



Furthermore, we caution that mass estimates derived from 
modelling the soft-excess are likely very highly uncertain. Gier- 
linski and Done (2004) in a systematic study of all radio-quiet PG 
quasars observed with XMM-Newton, found that the temperature 
of the soft excess was remarkably constant despite the 3 orders of 
magnitude range in black hole mass (their Figure 1, panel b) and 
order of magnitude range in L/L 0( jd. This apparent insensitivity of 
the soft-excess to either mass or luminosity has been noted before 
(Walter and Fink 1993; Czerny et al 2003; Porquet et al. 2004), 
and casts strong doubt on the interpretation of soft excesses as an 
accretion disc signature in QSOs. 

If the same uncertainty is applied to ULXs, we may require 
a different explanation for the origin of the soft excess. One pos- 
sibility is an optically thick, low temperature photosphere. This 
phenomenon should naturally arise in sources accreting at or close 
to the Eddington limit (King & Pounds 2003) without appealing 
to large black hole masses. In the King and Pounds model, any 
black hole accreting at close to Eddington will produce a strong 
Compton-thick outflow with an effective photospheric size of a few 
tens of .Rschw from the central source. Using their equation 18, and 
scaling down their numbers for the PG quasar 1211+143, then for 
a 1OM black hole accreting at close to Eddington we predict an 
outflow rate of ~ 10~ 7 M Q /yr, which implies an effective black- 
body temperature for the photosphere of a few x 100 eV, similar 
to the value we measure for Ho II X-l. However, we note that the 
largest challenge for this model is to explain how we still observe 
the hard X-ray component through the Compton-thick outflow (cf. 
Miller, Fabian & Miller 2004). 

If we return to the standard assumption - namely, the soft ex- 
cess does originate in the thermal emission of an accretion disc 
- then the mass estimates assume that the disc is in the "high" 
spectral state (i.e. radiating from the innermost stable circular or- 
bit around the black hole, = 6R g for a Schwarzschild black hole). 
We might then expect the rest of the X-ray spectrum and its vari- 
ability characteristics to conform to this state. Unlike some other 
ULXs, Ho II X-l does display the steep power-law continuum (in 
the empirical power-law + DISKBB fit), F > 2.4, characteristic 
of the high state (cf. Roberts et al. 2005). However, it does not 
display the correct variability characteristics. This is a particular 
problem as if an IMBH is present, even with the reduced mass de- 
rived above, the X-ray luminosity implies it is apparently accreting 
at ~ 0.1 — 0.5 L E dd- This accretion rate regime is observed in both 
AGN and XRBs, all of which show a substantial level of variability 
in the high state. Assuming that accretion physics works similarly 
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for IMBHs, this apparent failure of XRB properties to "scale-up" 
to an IMBH is a severe problem. 

There is a second way in which this ULX appears "abnormal" 
compared to Galactic XRBs in the high state. Physical modelling of 
the source spectrum reveals its putative coronal component is prob- 
ably both optically-thick and cool (r ~ 4, kT ~ 3 keV), in stark 
contrast to the optically-thin, hot (r ~ 1, kT ~ 50 - 100 keV) 
coronae found around Galactic black holes. The coronal X-ray 
emission also quantitatively dominates the observed X-ray emis- 
sion, contributing > 80% of the 0.3 - 10 keV flux, unlike Galactic 
black holes in the high state (cf. M c Clintock & Remillard 2005 6 ) 
In fact, an X-ray spectrum dominated by coronal emission is more 
indicative of the very high state. Interestingly enough, an optically- 
thick medium could provide a plausible means to explain the lack 
of any accretion disc features (in particular an Fe K line) in the X- 
ray spectrum. However, these arguments again imply that this ULX 
does not appear to be operating as a scaled-up high state XRB. 
Hence models assuming a ~ 1000 MqIMBH in the high state un- 
derlies Ho II X- 1 face considerable new challenges on the basis of 
our data. 

Is there an alternative explanation for the spectral and tempo- 
ral characteristics of this source? The Galactic microquasar GRS 
1915+105 might hold the key. In particular, Zhang et al. (2000) 
have demonstrated that its X-ray spectrum can be modelled by a 
three-layer atmospheric structure, broadly similar to the solar at- 
mosphere. These three layers are: (i) a cold and optically-thick disc 
(kT ~ 0.2 — 0.5 keV); (ii) a "warm layer" overlaying the cold disc 
with kT ~ 1 — 1.5 keV and r ~ 10; and (iii) a hot, optically- 
thin corona (kT ~ 100 keV, r ~ 1), though this last component 
is not always present 7 . This model has striking similarities to our 
physical model, in particular a cool disc seeding a warm, optically- 
thick scattering medium. Our lack of variability could then be ex- 
plained by the non-detection (within our limited spectral range) of 
the optically-thin coronal component and/or the source being in 
the x-class of GRS 1915+105 behaviour. Indeed, as the x-class of 
GRS 1915+105 appears typical of the very high state (Zdziarski et 
al. 2005), as is a corona-dominated X-ray spectrum, then the mass 
limit derived from band-limited variability considerations may be 
appropriate. Crucially, when combined with the non-dependence of 
the black hole mass on the disc temperature in this model, this al- 
lows us a far smaller black hole than conventionally assumed for 
this ULX « 100 Mq). This in turn implies that Ho II X-l must 
be accreting at (or perhaps above) the Eddington limit. We note 
that theoretical models for ULXs as stellar-mass black holes ac- 
creting from high mass donor stars certainly permit sufficient mass 
transfer to fuel such high luminosities (Rappaport, Podsiadlowski 
& Pfahl 2005), and speculate that this excessive mass transfer may 
also result in enough material transfering to the corona to make it 
optically-thick. 

Therefore we conclude that it remains plausible that Ho II X-l 
may possess a black hole no larger than a few tens of solar masses 
radiating at or above the Eddington limit. 



6 A caveat must be that if an IMBH is present, then to provide a fair com- 
parison one would have to use a softer X-ray band than the 2-20 keV band 
these authors use to define their black hole emission states. One might argue 
that our 0.3 - 10 keV data provides this, and is still very corona-dominated. 

7 Broad-band spectral studies of GRS 1915+105 can require models far 
more complex than this simple model, see for example Done, Wardzinski 
& Gierlinski (2004); Zdziarski et al. (2005). 



4 SUMMARY 

We present an unprecedented deep observation of one of the best 
IMBH candidate ULXs, Ho II X-l. This data provides a chal- 
lenge to previous results suggesting the presence of a ~ 1000 
M IMBH, most notably in the lack of significant variability over 
the timescales probed by our data, and the possible detection of 
an optically-thick corona, both of which are inconsistent with the 
simple argument for IMBHs based on the "scaling-up" of a Galac- 
tic XRB. A plausible alternative is that Ho II X-l is behaving like 
GRS 1915+105 in its x-class very high state; in this case the lack 
of observed variability places a limit on the mass of the compact 
accretor in Ho II X-l at < 100 M . 



ACKNOWLEDGMENTS 

We would like to thank the anonymous referee for a careful reading 
of the manuscript and for a number of helpful suggestions which 
led to an improvement in the work presented here. 

M.R. Goad and TP Roberts acknowledge financial support 
from PPARC. P. Uttley is supported by an NRC Research Associ- 
ateship. This work is based on observations with XMM-Newton, an 
ESA science mission with instruments and contributions directly 
funded by ESA member states and the USA (NASA). 



REFERENCES 

Begelman M.C., 2002, Apl, 568, L97 

Belloni T, Klein-Wolt M., Mendez M., van der Klis M., van 

Paradijs J., 2000, A&A, 355, 271 
Colbert E.J.M., Ptak, A.F., 2002, ApJS, 143, 25 
Cropper M., Soria R., Mushotzky R.F., Wu K., Markwardt C.B., 

Pakull M., 2004, MNRAS, 349, 39 
Czerny, B., Nikolajuk, M. Rozahska, A. Dumont, A.-M. 2003, 

A&A, 412, 317. 

Dewangan G., Miyaji T, Griffiths R.E., Lehmann I., 2004, ApJ, 
608, L57 

Dickey J.M., Lockman F.J., 1990, ARA&A, 28, 215 

Done C, Wardzinski G., Gierlinski M„ 2004, MNRAS, 349, 393 

Ebisawa K., Zycki P., Kubota A., Mizuno T, Wataria K., 2003, 

ApJ, 597, 780 
Ebisuzaki T, et al., 2001, ApJ, 562, L19 
Edelson R., Nandra, K„ 1999, ApJ, 514, 682. 
Fabbiano G., Trinchieri G, 1987, ApJ, 315, 46. 
Ferrarese L., Pogge R.W., Peterson B.M., Merritt D., Wandel A., 

Joseph C.L., 2001, ApJ, 555, L79 
Gebhardt K., et al., 2000, ApJL, 543, L5 

Gierlinski M., Maciolek-Niedzwiecki A., Ebisawa K., 2001, MN- 
RAS, 325, 1253 

Gierlinski, M. and Done, C. 2004, MNRAS 349, L7. 

Hidalgo-Gamez A.M., Sanchez-Salcedo F.J., Olofsson K., 2003, 
A&A, 399, 63 

Kaaret P., Ward M.J., Zezas A., 2004, MNRAS, 351, L83 

King A.R., Davies M.B., Ward M.J., Fabbiano G., Elvis, M., 2001, 

ApJL, 552, L109 
King, A.R. and Pounds, K.A. 2003, MNRAS 345, 657. 
Kording E., Falcke H., Markoff S., 2002, A&A, 382, L13 
Lehmann I., et al., 2005, A&A, 431, 847 
Madau P., Rees M.J., 2001, ApJL, 551, L27 
Makishima K. et al., 2000, ApJ, 535, 632 
Markowitz A., et al„ 2003, ApJ, 593, 96 



8 M.R. Goad, T.P. Roberts, J.N. Reeves, P. Uttley 



M c Clintock J. E., Remillard R. A., 2005, in Compact Stellar X- 
ray Sources, ed. Lewin W. H. G., van der Klis M., Cambridge 
University Press (Cambridge), in press l astro-ph/0306213 I 

M c Hardy I. M., Papadakis I. E., Uttley P., Page M., Mason K., 

2004, 348, 783 

M c Hardy I. M., Gunn K. E, Uttley P., Goad M. R„ 2005, MN- 
RAS, 359, 1469 

Miller J.M., Fabbiano G, Miller M.C., Fabian A.C., 2003, ApJL, 
585, L37 

Miller J.M., Fabian A.C., Miller M.C., 2004, ApJ, 607, 931 

Miller N.A., Mushotzky R.E, Neff S.G., 2005, ApJ, 623, L109 

Mitsuda K., et al., 1984, PASJ, 36, 741 

Miyaji T., Lehmann I., Hasinger G., 2001, AJ, 121, 3041 

Morrison R., McCammon D., 1983, ApJ, 270, 1 19 

Pakull M.W., Mirioni, L., 2002, preprint 

<astro-ph/0202488> 
Pourquet, D., Reeves, J.N., O'Brien, P., and Brinkmann, W. 2004 

A&A 422, 85. 
Pottschmidt K., et al., 2003, A&A, 407, 1039 
Rappaport S.A., Podsiadlowski Ph., Pfahl E., 2005, MNRAS, 356, 

401 

Remillard R. A., Muno M. P., McClintock J. E., Orosz J. A., 2002, 
in New Views on Microquasars, ed. Ph. Durouchoux, Y. Fuchs, 
and J. Rodriguez, Center for Space Physics: Kolkata (India), 49 
I astro-ph/0208402 1 

Roberts T.P, Warwick R.S., 2000, MNRAS, 315, 98. 

Roberts T.P., Warwick R.S., Ward M.J., Goad M.R., Jenkins L.P, 

2005, MNRAS, 357, 1363 

Soria R., Cropper M., Pakull M., Mushotzky R., Wu K., 2005, 

MNRAS, 356, 12 
Soria R., Motch C, Read A.M., Stevens I.R., 2004, A&A, 423, 

955 

Stobbart A., Roberts T.P, Warwick R.S., 2004, MNRAS, 351, 
1063 

Titarchuk L., 1994, ApJ, 434, 570 
Uttley P., M c Hardy I. M., 2005, in prep. 

Uttley P., McHardy I.M., Papadakis I.E., 2002, MNRAS, 332, 
231. 

Vaughan S., Iwasawa K., Fabian A.C., Hayashida K., 2005, MN- 
RAS, 356, 524 

van der Klis M., 2005, in Compact stellar X-ray sources, ed., 
W. H. G. Lewin, M. van der Klis, Cambridge University Press 
(Cambridge), in press l astro-ph/0410551 1 

Walter, R. and Fink, H.H. 1993, A&A 274, 105. 

Wilms J., Allen A., McCray R., 2000, ApJ, 542, 914 

Zdziarski A.A., Gierlinski M., Rao A.R., Vadawale S.V., Mikola- 
jewska J., 2005, MNRAS, 360, 825 

Zezas A.L., Georgantopoulos I., Ward M.J., 1999, MNRAS, 308, 
302 

Zhang S.N., Cui W., Chen W., Yao Y, Zhang X., Sun X., Wu X.- 
B., Xu H., 2000, Science, 287, 1239 



